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LATERAL STABILITY CHARACTERISTICS BETWEEN 
MACH NUMBERS OF 0.80 AND 1.57 AND SIMULATION OF COUPLED 
MOTION AT MACH NUMBER 1.30 OF A ROCKET -PROPELLED MODEL 
OF AN AIRPLANE CONFIGURATION HAVING THIN 
HIGHLY TAPERED 45° SWEPTBACK SURFACES 
By Charles T. D'Aiutolo and Allen B. Henning 

SUMMARY 


A rocket-propelled model of an airplane configuration having thin, 
highly tapered, 45° sweptback surfaces was tested in the Mach number 
range of 0.80 to 1.57 to obtain the lateral stability characteristics of 
the configuration. 

The results of the analysis of the Dutch roll oscillations indi- 
cated that from Mach number 1.3 to 1.57, the lateral-force derivative 
decreased and the directional stability was relatively constant. At a 
Mach number of 1-50, the value of the directional- stability derivative 
as obtained from the application of the time-vector method agreed with 
the value as obtained from the single-degree-of-freedom method. 

At a Mach number of 1.3 (a dynamic pressure of 2,060 pounds per 
square foot) , the model experienced a violent combined lateral- 
longitudinal motion, and large angles of attack and sideslip together 
with large rolling velocities were recorded. Approximate simulation 
of this combined motion on an analog computer using nonlinear equation 
of motion with five degrees of freedom was made and indicated possible 
values of the stability derivatives of the model while performing the 
roll-induced excursions in angle of attack and angle of sideslip. 


INTRODUCTION 


A systematic investigation of the dynamic lateral stability charac- 
teristics of airplane configurations at high subsonic, transonic, and 
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supersonic speeds is being conducted by the Langley Pilotless Aircraft 
Research Division by means of rocket-propelled models. Initial results 
of this investigation have been published in references 1 and 2 where 
the dynamic lateral stability characteristics of a swept-wing configu- 
ration (ref. 1) and a straight-wing configuration (ref. 2) are presented. 

This report contains the results of a flight test to determine the 
dynamic lateral stability characteristics of a rocket-propelled model 
of an airplane configuration having thin, highly tapered, k^>° sweptback 
surfaces. The Mach number range covered in this test was from 0.80 to 
1.57 nnd corresponds to a Reynolds number range (based on wing mean 
aerodynamic chord) of 6.7 X 1C)6 to 13-5 X 10^, respectively. The model 
was flown at the Langley Pilotless Aircraft Research Station at Wallops 
Island, Va. 

Stability derivatives were determined by analyzing the oscillations 
in yaw by means of the time-vector method as applied to rocket-propelled 
models. (See refs. 1 and 2). These oscillations were produced by the 
firing of pulse rocket motors. 

During the test, the model experienced a violent lateral -longitudinal 
coupled motion and large angles of attack and sideslip were recorded at 
high roll rates. A simulation study of this coupled motion was perforned 
on an analog computer to determine the stability derivatives of the model 
during this maneuver. 


SYMBOLS AND COEFFICIENTS 


In this paper the forces and moments acting on the model are 
referred to the body system of axes shown in figure 1 and may be trans- 
posed to the stability system of axes by use of the expressions derived 
in reference 3 • The origin of the axes system was at the center of 
gravity of the model which was coincident with the 0.193 mean aerody- 
namic chord of the wing. The symbols and coefficients are defined as 

total damping factor (logarithmic decrement of Dutch roll 
oscillation defined as being a positive number for a damped 
oscillation) 

acceleration along Y-reference axis as obtained from accel- 
erometer, positive to the right. 

wing span, ft 

horizontal-tail span, ft 


follows : 
a 


Ar/g 

b 
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mean aerodynamic chord, ft 

moment of inertia about body X-axis, slug -ft 2 

moment of inertia about body Y-axis, slug-ft 2 

moment of inertia about body Z-axis, slug-ft 2 

product of inertia in XZ-plane referred to body axes system 
(positive when the positive direction of the X-princlpal 

axis is inclined below the reference axis.) 
rolling moment, ft-lb 
pitching moment, ft-lb 
yawing moment, ft-lb 
Mach number 
mass of model, slugs 

period of the Dutch roll oscillation, sec 

angular velocities about X, Y, and Z body axes, radians/sec 
dynamic pressure, ipV 2 , lb/sq. ft 

Reynolds number 
total wing, area, sq. ft 
velocity, ft/sec 
model weight, lb 


nondimens ional lateral coordinate with respect to the tail 
span 

nondimens ional vertical coordinate with respect to the wing 
span 


air density, slugs/cu ft 
relative density factor, 


m 

pSb 


NA.CA RM L56A17 


natural circular frequency of the Dutch roll oscillation, 

radians/sec 

angle of attack, deg or radians 

angle of sideslip, deg or radians 

downwash angle, deg 

angle of pitch, deg 

angle of roll, radians 

angle of yaw, radians 

angle of sidewash, deg 

phase angle, deg 

sweep angle of quarter -chord line, deg 
incremental change in quantity 

square of the nondimens ional undamped natural frequency in 


yaw, 


Iz/ 


square of the nondimensional undamped natural frequency in 


pitch. 


c m a ( l sg 

I Y P 2 


fraction of critical damping ratio in yaw 


fraction of critical damping ratio in pitch 
lift coefficient, Lift/q^ 

lateral-force coefficient, Lateral force/q^S 
pitching-moment coefficient, Pitching moment/q^Scy- 
yawing -moment coefficient. Yawing moment/qoSb 
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C 2 

C np 

C Yp 

C *3 


Cv 

1 P 


c Y r 


C 


2 r 



rolling -moment coefficient. Rolling moment /c^Sb 
directional-stability derivative, dC n/ /dp, per radian 
lateral-force derivative, dCy /d(3, per radian 
effective dihedral derivative, per radian 


rate of change of yawing-moment coefficient with rolling- 

dC n 

angular-velocity factor, — ■ ■■-, per radian 


rate of change of lateral-force coefficient with rolling - 

SCy 


angular-velocity factor. 


damping-in-roll derivative, 


bC-, 


, per radian 


, per radian 


rate of change of yawing-moment coefficient with yawing- 

c>C n 

angular-velocity factor, j- . , per radian 

7 


rate of change of lateral-force coefficient with yawing- 

dCy 

angular-velocity factor, . * ■, per radian 

s (i?) 


rate of change of rolling-moment coefficient with yawing- 

dC 7 

angular-velocity factor, — , per radian 


<i) 


rate of change of yawing-moment coefficient with rate of 

ac n 

change of angle-of-sideslip factor, . v , per radian 

af 
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C Yr 


rate of change of lateral-force coefficient with rate of 


change of angle-of-sideslip factor, 


, per radian 




rate of change of rolling-moment coefficient with rate of 


change of angle-of-sideslip factor. 


dC- 


\2vy 


-, per radian 


"' m p 


'La 


rate of change of pitching-moment coefficient with angle of 
sideslip, dCm/dp, per radian 

lift-curve slope, dCL/da, per radian 


-'fflq 




rate of change of pitching-moment coefficient with pitching- 


bC. 


m 


, per radian 


angular-velocity factor. 


slope of pitching-moment curve, dC m /cfcx, per radian 




rate of change of pitching-moment coefficient with rate of 
change of angle-of-attack factor — r, per radian 

m 


drag coefficient, Drag/qoS 


Subscripts : 


w 

wing 

VT 

vertical tail 

PR 

pulse rocket 


The symbol I I represents the absolute magnitude of the quantity 
and is always taken to be positive. A dot over a variable indicates the 
first derivative of the variable with respect to time. Two dots indicate 
the second derivative with respect to time. The second subscript symbol 
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of the phase angles is used as a reference. A positive sign associated 
with the phase angle indicates that the first subscript symbol leads 
the reference, whereas a negative sign indicates that the first sub- 
script symbol lags the reference. 


MODEL, INSTRUMENTATION, AND TEST DESCRIPTION 
Model Description 


The general arrangement of the model is shown in figure 2, and the 
geometric and mass characteristics of the model axe given in table I. 
Photographs of the model and the model-booster combination are shown in 
figure 3* 

The fuselage was essentially a body of revolution and consisted of 
an ogival nose section, a cylindrical body section with rather large 
wing-root fairings for structural purposes, and a straight tapered 
afterbody section. The nose section contained the telemeter; the center 
section contained the wing and necessary wing-root fairings; the after- 
body section contained the horizontal and vertical tail surfaces. The 
fuselage was constructed of aluminum alloy with magnesium skin. 

The wing of the model was mounted along the fuselage reference 
line. It was made of aluminum and had NACA 65AOOL airfoil sections 
parallel to the free stream. It incorporated ^5° of sweepback along 
the quarter- chord line, and had an aspect ratio of 4 and a taper ratio 
of 0.2. The geometry of the horizontal tail was identical to the wing 
and it was mounted in the wing- chord plane extended. 

The vertical tail had an aspect ratio of 1.59 (based on the area 
to the fuselage center line), and contained NACA 65A004 airfoil sections 
parallel to free stream. It incorporated L5 0 of sweepback along the 
quarter-chord line, and had a taper ratio of 0.17* 


Instrumentation 

The model contained a standard NACA ten-channel telemeter. Measure- 
ments were made of the normal, longitudinal, and transverse accelerations 
near the center of gravity of the model. Rolling velocity was measured 
by means of a gyro-type instrument alined so that it was sensitive to 
velocities about the X-reference axis. The angle of attack and angle 
of sideslip were measured by an air-flow direction indicator located on 
a sting forward of the nose of the model. An air-flow direction indi- 
cator was also located in the plane of symmetry, above and ahead of the 
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horizontal tail and slightly forward of the vertical tail so that the 
flow angularities at a point near the vertical tail could be measured. 
(See figs. 2 and 3-) Total pressure was measured by a tube located on 
a strut mounted on the underside of the cylindrical section of the fuse- 
lage and by a tube located at a point above and to one side of the 
tapered afterbody section of the fuselage between the horizontal and 
vertical tails (see figs. 2 and 3) • 

In addition to measurements by the telemeter system, additional 
ground instrumentation was used. A rollsonde was used to measure the 
rate of roll of the model by means of polarized radio waves and the 
position of the model in space was determined by use of an MCA modified 
SCR 584 tracking radar set and the velocity of the model was obtained 
by use of the CW Doppler velocimeter radar set. Atmospheric data were 
obtained from a radiosonde released immediately before the model flight. 
Fixed and tracking motion-picture cameras were used to observe the 
condition of the model during most of the flight. 


Preflight Tests 

Prior to flight testing the mass and inertia characteristics of 
the model were determined (see table I) . The model was also suspended 
by shock chords and shaken by means of an electromagnetic shaker. This 
preflight test was performed to determine the structural natural fre- 
quencies of the model and these characteristics are also shown in 
table I. 


Flight Test 

The model was boosted to a Mach number of 1.62 and upon burnout of 
the booster rocket motor the model separated from the booster. After 
complete separation at a Mach number of 1.57 the model was disturbed 
in yaw by a series of four pulse-rocket motors. These rocket motors 
were located in the tapered afterbody section of the model (see fig. 2) 
and provided thrust normal to the XZ-plane. The pulse rocket motors 
which had about 20 lb-sec impulse with a normal burning time of 
0.06 second were timed to fire in sequence during the decelerating 
portion of the flight. The time histories of the resulting model motions 
were obtained by means of the MCA telemeter and instrument system. 

The flight conditions of the model are presented in figure 4 where 
the variation of air density, velocity, dynamic pressure, and relative 
density factor with Mach number are shown. These quantities are pre- 
sented so that a possible correlation of the data obtained from this 
test with data obtained from other tests may be made. The range of 
the Reynolds numbers of the present test is shown in figure 5» 


NACA RM L56AI7 


9 


ACCURACY AMD CORRECTIONS 
Accuracy 


The estimated probable errors in the measured quantities are shown 
in table II. Some or all of these measured quantities are reflected 
directly into the lateral stability derivatives Cy^, C n p, Czj3> C^p, 


and C nr 



The probable error in any one of the above derivatives 


due to all of the probable errors in table II was determined by the method 
shown in reference 2 , and at a Mach number of 1.50 the accuracies of 
the derivatives are: 3 percent for Cy^, 7 percent for C n p and C 


13 percent for C^, and 25 percent for C nr - C n ^ . The absolute values 


of these derivatives are as good as or better than the values indicated 
above . 


Corrections 

The readings of the air-flow direction indicator mounted in front 
of the model were corrected to the model center of gravity in order to 
determine the angles of attack and sideslip. Also the readings of the 
air-flow direction indicator mounted slightly forward of the vertical 
tail were corrected for model rotations so that when these readings were 
compared with the readings of the air-flow direction indicator mounted 
in front of the model (corrected to the center of gravity), the flow 
angularity at a point near the vertical tail could be determined. These 
corrections were made by considering flight-path curvature effects as 
shown in reference 4. 

The accelerometers were mounted as close as possible to the center 
of gravity of the model and values obtained from these instruments were 
corrected for the model motions in order to determine the accelerations 
of the center of gravity of the model. 

Frequency-response corrections to all instruments were not neces- 
sary and the model natural frequency was less than 4 percent of any 
instrument natural frequency . 
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RESULTS AND DISCUSSION 
Time History 


The time history of the lateral motions experienced by the model 
due to the yaw disturbances produced by the firing of the pulse rocket 
motors is shown in figure 6 where the variation with time of the Mach 
number, rolling velocity, lateral force, angle of sideslip, and the 
angle of attack are presented. 

The model separated from the booster at 3*45 seconds corresponding 
to a Mach number of about 1.57> which resulted in an angle of attack as 
well as a yaw and a roll disturbance to the model. (See fig. 6.) It 
is believed that this angle- of- attack change was small enough to have 
a negligible effect on the lateral motions. At t = 4.14 seconds 
(M = 1.44) the first pulse rocket fired and an induced angle- of-attack 
oscillation was recorded. This angle- of- attack oscillation quickly 
subsided and the model oscillated essentially in yaw. The gyro-type 
instrument used to measure rolling velocity became inoperative immedi- 
ately upon the firing of the first pulse rocket and accurate values of 
the rolling velocity could not be determined after this time. The 
second pulse rocket motor fired at t = 5*12 seconds (M = 1.30) and a 
violent combined lateral-longitudinal motion occurred in which an angle 
of sideslip of about 8° was recorded and an angle of attack greater 
than the 7° instrument range was attained. In order to determine the 
rate of roll during this coupled motion, the readings of the rollsonde 
were used. (See ref. 5 .) At best, the rollsonde can determine the 
magnitude of the roll rate, but does not accurately determine the phase 
between the roll rate and the angle of sideslip. A comparison was made 
of the roll rate as measured by the gyro-type instrument before it 
failed and the rollsonde to determine the accuracy of the magnitude of 
the roll rate as determined by the rollsonde. This comparison is shown 
in figure 6 (note that the gyro rolling- velocity scale is magnified). 
Reference 5 indicates that the rate of roll as determined from a 
rollsonde is least accurate when the rate of roll is changing rapidly. 
On this basis then, the rate of roll data shown in figure 6 when the 
rate of roll is changing rapidly is subject to considerable error and 
should be interpreted as showing only a trend. Therefore, during the 
combined lateral-longitudinal motion the roll rate built up rapidly and 
a value of about -45 radians/second was attained. 

At t = 5.40 seconds (M = 1.24) the third pulse rocket fired and 
the angle of sideslip decreased rapidly and oscillated at relatively 
low amplitude. The rate of roll increased somewhat and then decreased 
considerably, but the angle of attack remained above the limit of the 
instrument range (7°). (See fig. 6.) The model was still experiencing 
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a combined lateral-longitudinal motion as shown by the frequency of the 
lateral oscillations and the magnitudes of the angle of attack and 
rolling velocity. 

The fourth pulse rocket fired at t = 6.67 seconds (M = 0.88) and 
the rate of roll decreased rapidly. Also the angle of attack decreased 
rapidly and tended to subside . At this time the model experienced a 
lateral disturbance with very little pitch coupling. 

From the time history of the lateral motions without pitch coupling 
present, the periods of the Dutch roll oscillation were determined and 
are shown in figure 7* These periods were determined from the oscilla- 
tions in angle of sideslip, lateral force coefficient, and roll velocity. 
In the Mach number range where pitch coupling was present, the periods 
were calculated by two methods . The first method considered only the 
basic frequency, whereas the second method considered the rolling fre- 
quency as well as the basic frequency. (See ref. 6.) These data are 
also shown in figure 7 for Mach numbers at which the roll frequency 
could be accurately determined. The trend of the data indicates that 
the period decreases with increasing Mach number. 


Trim Characteristics 

The variations of the model trim characteristics with Mach number 
are shown in figure 8 as trim angle of attack, trim angle of sideslip, 
and trim lateral-force coefficient. These characteristics were determined 
when pitch coupling was negligible. The trim angle of attack remains 
relatively constant at subsonic and supersonic speeds at a value of 0.80°. 
The trim angle of sideslip is negative at subsonic speeds and becomes 
positive at supersonic speeds, while the trim lateral-force coefficient 
remains negative at both subsonic and supersonic speeds. 


Lateral-Force Derivative 

Two typical cross plots showing the variation of lateral-force 
coefficient with angle of sideslip are presented in figure 9* These 
cross plots are shown for M = 0.86 and M = l.k2 and indicate little 
difference in the slope of the lateral-force coefficient against angle 
of sideslip, Cyp. Presented in figure 10 is the variation of Cy^ with 

Mach number which was determined from the cross plots of Cy against 3. 
The variation of Cy^ with Mach number is presented only in the Mach 

number range when pitch coupling was negligible. Also shown in figure 10 
is the vertical- and horizontal-tail- contribution to Cyp calculated from 

the theory presented in reference 7. These results, along with an 
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estimate of Cy^ of the wing-fuselage of about -0.1 per radian, indi- 
cate that at supersonic speeds the main contribution to Cy^ is the 

vertical tail. The level of the measured Cy^ is in good agreement 

with the calciilations when due consideration, on a qualitative basis 
only, is given to possible effects such as aeroelasticity . 


Vector Plots 


As was previously mentioned, the gyro-type instrument used to 
measure rolling velocity became inoperative early in the flight of the 
model. During the major portion of the flight the model experienced a 
violent combined lateral -longitudinal motion; therefore, it was only 
possible to apply the vector method of analysis to the lateral motions 
at one Mach number, M = I. 50 . These vector plots of the lateral equa- 
tions of motion are presented in figure 11 and were obtained by the 
method presented in reference 1. The analysis was performed for an 
undamped natural circular frequency, 00, of 38 . 16 radians/sec; a total 
a, of 1.8l per sec; 

JL 

ipr 


damping factor, 
velocity to angle of sideslip, 


an amplitude ratio of rolling 
radians/ sec t 
degree ’ 


of 157.22 


and the 


phase of the rolling velocity to angle of sideslip, ® 0 f 131.56°. 


Shown in figure ll(a) is the vector diagram of the lateral force 
equation. The solution of this equation yielded values of -|i| and 

which were used in the vector diagrams of the rolling and yawing 

moment equations. Figure ll(b) is the vector diagram of the rolling 
moment equation. In order to solve this equation, it was necessary to 

assume values for either Cj r or . Values of Cj r were assumed 

P 

and values of Cjp and Cj were determined as functions of C^ r : 

C ir = 0.10, 0.15, and 0.20. Shown on figure 11(b) are the vector 

lengths Czp-M, %(|jr)!fp C ^r(^)jtf* Listed beloV are the 

the values of and obtained from these vectors. 
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For values of 


c Zr = 0.10 

c Zp = -0.078 

c Zp = -0.251 

Cz r = 0.15 

C Zp = -0.076 

c lp = -0.256 

C Zr = 0.20 

c Z p = -0.071*- 

C 7 = -0.221 

ip 


The vector diagram for the yawing moment equation is presented in fig- 
ure 11 (c) and values of C n p and ^C nr - C n ^ were determined as func- 
tions of C np ; C np = -0.05, 0, and 0.05. Shown on figure ll(c) are the 
vector lengths C np -j||, (C nr - -jjj-, and C np ^|||. Listed 

below are the values of C n p and obtained from these vectors. 

For values of 


Cn p 

= -0.05 

C np =0.171 

(°n r " c np) 

= -0.510 

Cn p 

= 0 

C np =0.178 

(C nr " C np) 

= -0.151 

Cn P 

= 0.05 

Cnp =0.185 

( C n r “ C np) 

= 0.015 


A comparison of C np as determined from the vector method and the 
single -degree -of -freedom method is presented in figure 12. The single- 


degree-of-freedom value for C 


np 


is determined directly from the 

kit 2 I<7 

periods of the Dutch roll oscillations, C nR = ^ 

P qoSbP" 


The agreement 


between C n p, as determined from the single-degree-of -freedom method 

and the vector method, when C np = 0 is good. The single-degree-of- 

freedom method for the determination of C n p shows that C n p is nearly 

constant at supersonic speeds. Also shown in figure 12 is the vertical- 
and horizontal-tail contribution to C n p calculated from the theory 

presented in figure 8 . Available data on wing -body combinations similar 
to that of the model reported herein (except for the large wing-root 
fairings) indicate values of C n p of about -0.1 per radian. Agreement 


between C n p 


for the complete configuration based on these data and the 


measured data is good, considering possible effects such as loss in 
effectiveness of the vertical tail due to aeroelasticity . These results 
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are in agreement with the results discussed previously with regard to 
the lateral force derivative, Cy^ • 

Simulation of Combined Lateral-Longitudinal Motion 

Si mill at ion . - A simulation of the combined lateral-longitudinal 
motion experienced by the model at M = 1.3 and q^, = 2,060 lb/sq. ft 

was attempted in order to determine the stability characteristics of 
the model while undergoing this maneuver. The model was assumed to have 
five degrees of freedom and the forces and moments acting on the model 
are described by the following equations of motion: 

Rolling : 


P = 


% - + (iJBU + fiXZApq + u, p 


L x 


^lc 7P + 

2VI X 


c + C lR (a)p + ^ c 


2VI X 


IX P ix 


'ipR 


Pitching: 


<1 = 


' T Z - Ix \pr + (5$ r 2 - f^Av 2 + C m „d + 




2VI y 


Yawing: 


2o£! Cnw i t i£! + S & Cne p 

2VIy » Iy «- Iy p 




^^ ■-( c n r " Cn p) r + c np( a )P + c n 


2VI Z 


iz 
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Y- force : 


0 = -r + ap + Cy r 0 + Cy^ 
mV P mV PR 


Z-force: 


d = q. - pp - ^/c L + Cp\a 
mV \ a / 

The terms C^(a) and C n p(a) denote that C^p and C n ^ are variant 
with angle of attack. The terms Cjp R , C n _^, an(i ^YpR respec- 

tively, the rolling-moment coefficient, yawing-moment coefficient, and 
the lateral-force coefficient produced hy the firing of the small pulse 
rocket motors . 


Certain simplifying assumptions were made in the above equations. 
The gravity terms were neglected in the Y- and Z-force equations since 
the flight conditions of the model resulted in values of these terms 
which were very small. Aerodynamic terms in the Y-force and rolling 
moment equations due to rate of change of angle of sideslip were assumed 
to be zero, and in the yawing moment equation it was assumed ^ = (3 so 
that C n ^ and C n £ may be combined into one term. Further, it was 

assumed that in the Y-force equation the only aerodynamic term retained 
was that due to angle of sideslip since this quantity as measured con- 
tained the Y-force contributions due to ^ and $ as well as (3 . 

Since these equations of motion are nonlinear, an analog computer 
was employed in order to determine a solution. The aerodynamic char- 
acteristics of the model were estimated from various rocket model tests, 
wind-tunnel tests, and theory, whereas the mass characteristics used 
were those measured in the preflight tests of the model. The procedure 
employed in the simulation was to adjust systematically the estimated 
aerodynamic characteristics of the model until the solution of the equa- 
tions of motion matched the flight test record of the combined motion. 

A considerable number of runs were made on the analog computer and 
changes in values of C^, C n ^ (from estimated negative value of C n ^ 

to positive value of C n p), C^, or had- large effects on the 

simulation, while wide variations in Cyp, C^, C^, and had 

little effect. The best simulation obtained appears in figure 13 . Zero 
time in this figure corresponds to the flight time when the second pulse 
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rocket was fired (t = 5*12 sec, see fig. 6) . In order to achieve this 
simulation of the flight test record, a linear* variation of and 

Cnp with angle of attack had to be incorporated into the equations of 
motion, but all other derivatives were considered invariant with angle 
of attack. The stability derivatives used are shown in table III. 

The simulation shown in figure 13 is qualitatively good. Quanti- 
tative differences between the flight test results and the analog results 
indicate the need for more detailed studies of this type before any 
conclusions can be reached as to the accuracy of the derivatives obtained 
by this method. 

Steady rolling analysis .- An analysis was made (based on ref. 6) 
of the combined lateral-longitudinal motion in order to determine how 
effective the method is in predicting the roll rate required for diver- 
gence. Values of the stability derivatives determined from the simula- 
tion were used in this ^analysis and the critical damping ratio in pitch 
was determined from a two-degree -of -freedom analysis and the critical 
damping ratio in yaw was also determined from a two-degree -of -freedom 
analysis. The results appear in figure Ik and indicate that the model 
would experience a divergence in angles of attack and sideslip when the 
rolling rate is between 39 and 4-2 radians/second. It was impossible to 
determine exactly the roll rate the model experienced when divergence 
occurred (see fig. 6), but the model did diverge at a roll rate consid- 
erably less than the maximum value recorded (-‘4*5 radians/sec). 

It appears then that the analysis of reference 6 does predict that 
a divergence would occur, but due to the inaccuracy of the measured roll 
rate it is impossible to tell just how accurate the analysis is in pre- 
dicting the exact value of the roll rate required for divergence. Also, 
as previously mentioned, the analysis accurately predicts the basic 
frequency of the model under steady rolling conditions . 


Downwash and Sidewash 

Data that have been obtained from an attempt to measure downwash 
and sidewash are presented herein to show what type of data can be 
expected from the measuring system used on this model. The method used 
to measure downwash and sidewash employed the nose air-flow indicator 
and an additional air-flow indicator that was sting mounted to the 
vertical tail. (See figs. 2 and 3.) The tail air-flow indicator meas- 
ured the flow direction at a point which was in the plane of symmetry, 
above and ahead of the horizontal- tail surface. 

The angular difference in the angle of attack and angle of sideslip 
between the free stream or nose air-flow indicator and the tail-mounted 
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air-flow indicator represents the angle of downwash and the angle of 
sidewash induced by the flow of air over the wings and fuselage . 

The downwash data which were obtained from, this test are shown in 
figure 15 and compared with various values obtained from previous experi- 
ments, theories, and calculations of references 8 to 12. The error in 
de/da was calculated by the method in reference 15 and was found to be 
from 12 to 14 percent of the measured value . The measured values of 
de/da from this test were taken from the largest amplitude of the angle 
of attack for each readable oscillation. The largest amplitude permitted 
the least amount of error in de/da. The values obtained from this test 
are somewhat high in comparison to the reference values . It should be 
pointed out that the values obtained from this test were values for one 
specific point in the downwash field and not the average value over the 
horizontal tail; therefore, this higher value of de/da could possibly 
be justified. 

The sidewash, da/dp, was determined by using the same procedure 
as for downwash. One test point was determined and its value was 0.19 
at a Mach number of 1.412. The error involved was on the order of 
26 percent of the measured value. This sidewash value is believed to 
be influenced entirely by sideslip and not by roll, as the roll parameter, 
pb/2V, was less than 0.003 at this Mach number. No comparison values 
for this type of sidewash were found; therefore, these data are included 
herein purely for academic reasons. 

The dynamic pressure at a point above the tapered afterbody section 
of the fuselage between the horizontal and vertical tails was determined 
from a comparison of the measured values of the free-stream total head 
pressure and the total head pressure measured at this point by consid- 
ering the static pressure to be constant and equal to the free-stream 
value. The data indicate that within the limits of accuracy of the 
measured quantities the dynamic pressure at a point above and to one 
side of the tapered afterbody section of the fuselage between the hori- 
zontal and vertical tails was equal to the free-stream dynamic pressure. 

(a < 1°; 3 < 1°.) 


CONCLUDING REMARKS 


From the results of a flight test to determine the lateral stability 
characteristics of a rocket-propelled model of an airplane configuration 
having thin, highly tapered, 45° sweptback surfaces, the following con- 
clusions are indicated: 

At a Mach number of 1.3 the model experienced a violent combined 
lateral-longitudinal motion and large angles of sideslip and attack 
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together with large rolling velocities were recorded. Approximate simu- 
lation of this combined motion on an analog computer using nonlinear 
equations of motion with five degrees of freedom was made and indicated 
possible values of stability derivatives of the model while undergoing 
this maneuver. 

The results from the analysis of the Dutch roll oscillations indi- 
cated that from a Mach number of 1.3 to a Mach number of I .56 the lateral- 
force coefficient decreased and the directional stability was relatively 
constant with Mach number. It was only possible to apply the time-vector 
method of analysis to the lateral motions at one Mach number: a Mach 

number of 1.50. The value obtained for the directional- stability deriv- 
ative agreed with the value obtained by the single-degree-of- freedom 
method . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 6 , 1955- 
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TABLE I.- MASS AND GEOMETRIC CHARACTERISTICS 

Wing: 

Total area, sq ft 4.4l 

Span, ft 4.20 

Aspect ratio 4 

Taper ratio 0.2 

Sweep of quarter-chord line, deg 45 

Mean aerodynamic chord, ft 1.208 

Dihedral, deg 0 

Incidence, deg 0 

Airfoil section, parallel to free stream 65A004 

Horizontal tail: 

Total area, sq ft 0.88l 

Span, ft 1*875 

Aspect ratio 4 

Taper ratio 0.2 

Sweep of quarter-chord line, deg 45 

Airfoil section, parallel to free stream 65A004 

Vertical tail: 

Area, sq ft (total included area to center line of model) . . O.885 

Span, ft I.I87 

Aspect ratio 1*59 

Taper ratio 0.172 

Sweep of quarter-chord line, deg 45 

Airfoil section, parallel to free stream 65A004 

Fuselage : 

Length, ft 6.33 

Fineness ratio 10.72 

General : 

Weight, lb 99*75 

I X , slug-ft2 0.55 

Iy, slug-ft 2 6.23 

I Z , Slug-ft 2 6.79 

1X2, slug-ft 2 O.0655 

Center-of -gravity location 0.193c 

Wing first bending frequency, cps 46.5 

Wing second bending frequency, cps 160 

Horizontal tail first bending frequency, cps 112 


ro 

ro 


TABLE II.- ESTIMATED ACCURACY OF BASIC QUANTITIES 


[values shown are positive or negative quantities^] 


Mach 

number 

Accuracy of - 

w, 

percent 

dz, 

percent 

percent 

x xz> 

percent 

M, 

percent 

Qo> 

percent 

<0, 

percent 

111 

K/gl 

V 

deg 

a, 

deg 

a, 

sec - "*" 

|3| 

percent 

|3 1 ' 

percent 

1.50 

0.5 

2.0 

2.0 

16.0 

1.0 

2.0 

3.0 

3.0 

2.0 

3.0 

0.5 

0.1 


t 
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TABLE III.- VALUES OF DERIVATIVES AS OBTAINED FROM 
SIMULATION OF COUPLED MOTION. M = 1.3 


[All values for zero angle of attack. All 
derivatives in radian measure]] 


c np = 0.25 

c n r “ c np - 

0 

C-i 

TD 

II 

1 

O 

b 

00 

c 2r = 0.10 

C Y(3 = -0.76 

Cl* " 3.3 

C np = 0.20 

Cn^ = -0.86 

% = -0.32 

Cmq = _1+ .T 

Cjnp = 0 . 17 

ClT£ = 3.0 


C D = 0.035 


Note: 

that is, C 1 

C = -0.16; 

at a = 5°, 


and C n p were assumed variant with angle of attack, 
p increased with angle of attack until at a = 7°> 
while C n p decreased with increasing angle of attack and 
"np 
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Figure 1.- Sketch showing the body-axes system. Each view presents a 
plane of the axes system as viewed along the third axis. Arrows 
indicate positive direction of forces, moments, and angles. 




Figure 2.- General arrangement of model. All dimensions are in inches 

unless otherwise noted. 
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(a) Side view. L- 82 I 4 J+I.I 



(b) Top view. L- 82[|.59 • 1 


Figure 3-- Photographs of model and model-booster combination. 
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(c) Model-booster combination on launcher. 

Figure 3 .- Concluded. L- 82832.1 
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Figure 4.- Variation of relative density factor, dynamic pressure, 
velocity, and air density with Mach number. 


NACA RM L56A17 


NACA RM L56A17 


29 



Figure 5 .- Variation of Reynolds number with Mach number. 
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Figure 6.- Time history of model flight. 
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Flagged symbols represent data corrected 
for rolling frequency. (Reference 6). 
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Figure 7>- Variation of period of lateral oscillations vith Mach number. 
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Figure 8.- Variation of the trim characteristics with Mach number. 
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P,deg 


Figure 9*- Typical variation of side force coefficient with angle of 

sideslip. 
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Figure 10.- Variation of the lateral-force derivative with Mach number. 
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(a) Lateral- force equation. 


Figure 11.- Vector plots for the calculation of the lateral stability 

derivatives at M = 1 . 5 . 
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(b) Rolling-moment equation. 
Figure 11.- Continued. 
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Figure 11.- Concluded. 
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Figure 12.- Variation of the directional- stability derivative with 

Mach number . 


NACA RM L56A17 


NACA RM L56A17 


39 


Flight test 

Simulation 



Figure 13.- Comparison of flight test and simulation results of combined 
lateral- longitudinal motion. M = 1.3; q. Q = 2,060 lb/sq. ft. 
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Figure 14.- Effect of steady rolling on the pitch and yaw stabilities. 

M = 1.3. 
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Figure 15.- Comparison of downwash with other experimental data. 
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